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Brane dark energy cosmologies, leading to various possible evolutions of our universe, are investi- 
gated. The discussion shows that while all these models can be made arbitrarily close to the standard 
ACDM cosmology at present, their future evolutions can be very different, even diverge with time 
in a number of ways. This includes asymptotic de-Sitter evolution, Little Rip with dissolution of 
bound structures, and various possible singularities, as the Big Rip, a sudden future singularity 
(Type II), and Type III and Type IV cases. Specifically, some interesting effects coming from the 
brane tension are investigated. It is shown, in particular, that the Little Rip occurs faster on the 
brane model than in usual FRW cosmology. And in the asymptotic de-Sitter regime the influence of 
the brane tension leads to a deviation of the value of the effective cosmological constant from that 
O , corresponding to ordinary dark energy. As a consequence, the value of the inertial force from the 

accelerating expansion can greatly exceed the corresponding inertial force in ordinary cosmological 
' models. 
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I. INTRODUCTION 



The discovered cosmic acceleration [HQ can be interpreted in terms of the so-called dark energy (for recent reviews, 
see [3, IH), which in a way or other exhibits quite strange properties. However, some attempts at keeping General 
Relativity (GR) and vacuum fluctuations (sort of a Casimir effect) are still being considered, see e.g. [__, and as 
well known, with the help of an ideal fluid GR can be rewritten, in an equivalent way, as some modified gravity. 

The equation of state (EoS) parameter too for dark energy is known to be negative: 



w D = Pd/pd < , (1) 

[ being here pd the dark energy density and prj the pressure. According to astronomical observations, dark energy 
currently accounts for some 73% of the total mass-energy of the universe (see, for instance, Ref. [7|).Even if these 
observations consensually favor the standard ACDM cosmology, being actually rigorous the uncertainties in the 

. determination of the EoS dark energy parameter w are still too large to define which of the three (quite different) 
cases w < — 1, w = — 1, and w > — 1 (remaining always in the vicinity of w = — 1, of course) is actually realized in 
our universe, the present observational value (within the framework of the phantom energy model with constant Wq) 
being w = -1.04±° ^ 9 [| Q. From the theoretical point of view it is interesting to consider models with evolving wq, 

. £h 'i 

If w < — 1 (the phantom dark energy case [l(|) all four energy conditions are violated. According to quantum field 
theory a field of this kind is unstable. But the possibility to have the effective EOS w e s < — 1 for dark energy in 
scalar-tensor gravity without gho sts and other instabilities was shown in [ll|. Phantom dark energy can lead to a 
Big Rip future singularity (Tol Il2| - fl5l ]. where the scale factor becomes infinite at a finite time in the future. Another 
possible scenario is a sudden (Type II) singularity [l6[ where the scale factor is finite at the Rip time. However, if 
w tends to —1 sufficiently fast, then no finite-time future singularity occurs [l7j - l20j ]. Also, within the framework of 
braneworld models the possibility to have w c h < — 1 without future singularities was originally pointed out in (2lj j . In 
this case the energy density increases with time or remains constant. If the energy density grows, the disintegration 
of any bound structure eventually occurs, in a way similar to the case of the Big Rip singularity. 

The future evolution of our universe depends on the choice of the EoS satisfied by dark energy: 

p = g(p), (2) 
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with g a function of the energy density. 

The aim of this article is to develop a general method for the construction of dark energy models on the brane which 
are compatible with observational data and which could provide us with various possible scenarios of universe evolution. 
This question is analyzed from the viewpoint of the dark energy EoS and, in particular, from its corresponding 
description in terms of a scalar field theory. To this end, Sect. II contains a brief overview of the EoS formalism in 
FRW cosmology. Sect. Ill deals with the explicit generalization of the EoS formalism to the brane. The classification 
of cosmological models with different future evolutions is there undertaken. Non-singular models are considered in 
Sect. IV. It is demonstrated there that some of these models can actually be made indistinguishable from the ordinary 
ACDM cosmology, thus remaining stable for very long time before the disintegration of bound structures takes place. 
In Sect. V cosmologies which develop future "soft" singularities are considered. The properties of these models are 
seen to be actually similar to those of the previous section. A detailed numerical comparison of the predictions of the 
models with the observational data leads to the conclusion that they are indeed indistinguishable from the ACDM 
model in describing the universe evolution up to the present epoch. Moreover, as already mentioned, they may remain 
stable for billions of years before reaching a future singularity. In this sense, any of these models represents an actually 
viable alternative to ACDM. A key qualitative result of our study is that current observations make it essentially 
impossible to determine whether or not the universe will end up in a future singularity. The paper finishes with a 
section devoted to discussions. 



II. DARK ENERGY IN FRW COSMOLOGY: BRIEF OVERVIEW 



For a spatially flat universe, with metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (3) 



the cosmological equations are 



where p and p are the total energy-density and pressure, respectively, a is the scale factor, ' = d/dt, and we use natural 
system units in which 8ttG = c = 1 . 

The EoS for dark energy is conveniently rewritten as 

Pd = -Pd - /(ad), (5) 

where J(pd) is a function of the energy-density. The case f(po) > corresponds to w < —1, while the case f(p) < 0, 
to w > — 1. For simplicity, we omit the subscript D. If dark energy dominates, one can neglect the contributions of 
dark matter and radiation. Then, from Eqs. Q, one gets the following link between time and J(pd)'- 

' '' - VP- (6) 



For quintessence the energy density decreases with time, x < xq, while for phantom energy x > xq- Let us consider 
possible variants of the behavior of © (compare with [la]). 

1. The integral © converges for p — > oo. Therefore, a finite-time singularity occurs: the energy density becomes 
infinite at a finite time, tf. The expression for the scale factor 

(2 [ x xdx\ 

a = a exp - / — - , (7) 



3 Jx /O) 

indicates that there are two possibilities, namely: 

(la) the scale factor diverges at finite time (a "Big Rip"); 

(lb) the scale factor reaches a finite value and a singularity (p —¥ oo) occurs. It is of type III, in the notation 
of Ref. [33|]. The key difference between (la) and (lb) is that the energy density in the second case grows so 
rapidly with time that the scale factor does never reach the infinite value. 



3 



2. The integral ([6]) diverges for p —> oo (these models are described in detail in [13, The energy density grows 
with time but not quickly enough for a Big Rip occurrence. According to the terminology in (l7j we have a 
so-called "Little Rip" : eventually a dissolution of all bound structures, one after the other, will happen in the 
future (the first and very simple example of Little Rip was given in |22|). 

These scenarios can be realized only in the case of having a phantom energy. The next two variants are possible 
both for a phantom energy and for quintessence. 

3. The integral © diverges when p — > pf. The dark energy density tends asymptotically to a constant value (an 
"effective cosmological constant" [HI). These models constitute natural alternatives to the ACDM model, which 
also leads to a non-singular cosmology. Nevertheless, even for the non-singular asymptotically de Sitter universe, 
the possibility of a dramatic rip which may lead to the disappearance of bound structures in the universe still 
remains. 

4. Another interesting case occurs if f(x) — > ±oo at x = xj, i.e., the dark energy pressure becomes infinite at finite 
energy density. The second derivative of the scale factor diverges while the first derivative remains finite. It is 
interesting to investigate the family of phantom energy models with such kind of finite-time future singularities. 

An equivalent description in terms of a scalar theory can be obtained using the equations: 

p = ±tf/2 + V(& 1 p = ±4> 2 /2-V(0), (8) 

where <f> is a scalar field with potential V(<fi). The minus sign in front of the kinetic term corresponds to the phantom 
energy case. For the scalar field and potential one derives, respectively, the following expressions 

9 t" x dx 

d>(x)=d> ±^= -==, V(x)=x 2 + f(x)/2. (9) 
V3 Jx vl/wl 

Combining Eqs. (j9|) one obtains the potential as a function of the scalar field. 



III. DARK ENERGY MODELS ON THE BRANE: THE EQUATION OF STATE FORMALISM 

We now consider a spacetime which is homogeneous and isotropic along three spatial dimensions, being our 4- 
dimensional universe an infinitesimally thin wall, with constant spatial curvature, embedded in 5-dimensional space- 
time [lil, [24j]. In the Gaussian normal coordinate system, for the brane which is located at y = one gets: 

ds 2 = -n 2 dt 2 + a 2 (t, y^dx'dx 3 + edy 2 , (10) 

where 7^ is the maximally 3-dimensional metric, and e = 1, if the extra dimension is space- like, and e = —1, if it is 
time-like. Let t be the proper time on the brane (y = 0), then n(t,0) = 1. Therefore, one gets the FRW metric on 
the brane 

ds 2 y=Q = -dt 2 + ^(t^^dx'dx 3 . (11) 
The 5-dimensional Einstein equations have the form 

Rab - t^QabR = X 2 Tab + kg AB , (12) 

where A is the bulk cosmological constant. The next step is to write the total energy momentum tensor Tab on the 
brane as 

T£ = Sp(y), (13) 

with Sg = di&g(—pi,,Pb,Pb,Pb, 0), where pb and pb are the total brane energy density and pressure, respectively. 

One can now calculate the components of the 5-dimcnsional Einstein tensor which solve Einstein's equations. One 
of the crucial issues here is to use appropriate junction conditions near y — 0. These reduce to the following two 
relations: 

dn x 2 . X 2 da X 2 ^ 

ndy ly=0+ 3 2 ady\ 0+ 6 
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After some calculations, one obtains 

36 6 nZ " 4 



This expression is valid on the brane only. Here H = a(t,0)/a(t,0) and C is an arbitrary integration constant. The 
energy conservation equation is correct, too, 



Pb + 3% b +p b ) =0. (16) 
a 



Now, let pb = p + A, where A is the brane tension. For a fine-tuned brane with A = eA 2 x 4 /6, we have the following 
equation (for k — 0) 

h 2 eAx 4 /OA, p \ , C 



a 2 6 3 V 2A 

In this paper we consider a single brane model which mimics general relativity (GR) at present but differs from GR 
at late times. We set 87rG = eerx 4 /6. As is easy to check, two cases arise: e = 1 and A > 0, or e = — 1 and A < 0. 
For simplicity, we also set C — (the term with C is usually called "dark radiation"). In fact, setting C^O does not 
lead to additional solutions on a radically new basis, in the framework of our approach. Eq. (|17[) can be simplified to 

P -(l + ^ 
a 2 3 V 2A 



One can see that Eq. (fl8|) for p « IAJ differs insignificantly from the FRW equation. The brane model with a 
positive tension has been discussed in (2514271 ] in the context of the unification of early and late time acceleration eras. 
The braneworld model with a negative tension and a time-like extra dimension can be regarded as being dual to the 
Randall-Sundrum model 28-30]. Note that, for this model, the Big Bang singularity is absent. And this fact does 
not depend upon whether or not matter violates the energy conditions [3l[ . This same scenario has also been used to 
construct cyclic models for the universe [32| . 

The EoS is taken under the form It is convenient to put 

f 2Asinh 2 2 A > 0, e = 1, , in . 
^ = \2|A|sin 2 !, A<0, e = -l. (19) 

The parameter rj varies from — oo to oo, for A > 0, and from — 7r to 7T, for A < 0. Then, one can write equations 
similar to Eqs. © and ([7J 



3 J m HlY 



do exp 



A rv ,^pan) A > 0, e = 1, 
3 Jvo /(»;) J ' ' 



<V) = < >,,, , < (21) 

«oexp(f A<0, e=-l. 

The situation is very close to the one for the convenient FRW cosmology. In the case of a positive tension the following 
possibilities can be realized: 

1. If the integral (I2U1) converges while ([7]) diverges, we have a Big Rip. It is interesting to note that the Big Rip 
on a brane considered in [34| occurs faster than in the FRW cosmology. To illustrate this remarkable fact, one 
can rewrite Eq. (|20[) in the explicit form: 



2 f x 

t{x) = 7^ / 7" 

V3 Jx (i 



dx 



For the simplest EoS with constant EoS parameter wo = — 1 — a 2 , the function is f(x) = a 2 x 2 . If p >> A, then 
the dark energy density grows with time substantially faster than in ordinary cosmology (A — > oo). 

2. If the integrals ([20)1 and (|2"TT) diverge when r] — > oo, then a Little Rip occurs. The acceleration of the universe 
increases with time definitely faster than in the FRW universe, owing to the brane tension (see the corresponding 
time equation for the case (1)). 
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3. Asymptotic de Sitter expansion is realized if / —> for 77 — > 77/, and the integral (|20|) diverges. 

4. There is a type III singularity if both integrals converge when 77 — > 00. 

5. If f(rj) —?• 00 for 77 — 77 j , the universe ends its existence in a sudden future singularity. 
The case of negative tension allows for the following interesting possibilities: 

1. Asymptotic de Sitter expansion, if f(j)) — ¥ for 7/ — > r/j. 

2. Asymptotic breakdown (i.e. the rate of expansion of universe tends to 0) will occur if f(j]) — > for 77 — > tt. 

3. A sudden future singularity, if f(jj) — > 00 when 7/ — > 77/. 

For the scalar field and potential we get, respectively, the following equations 



^WWJttof i<S (22) 

One should note that dark energy with an EoS such that f(p) ~ p 7 , with 7 < 2, leads to a Big Rip on the brane 
while, in the case of the conventional FRW universe, such dark energy leads to a Little Rip only. 

Some additional remark is in order. Any realistic cosmological scenario should take into account that dark energy is 
not a single component of the universal energy. The key test for dark energy models consist in matching the modulus 
p vs redshift z — a^/a — 1 relation obtained by the SNe project. As is well known, 

H(z) =const + 51gL>(z). (23) 

The relation for the luminosity distance Dl(z) as a function of the redshift, in the FRW cosmology (FC), is 

D[ C = + z) I" [On(l + zf + Q D p(z)/p ] - 1 ' 2 dz . (24) 

"0 Jo 

Here, fl m is the total fraction of matter density, flo the fraction of dark energy energy density, and Hq is the current 
Hubble parameter. For simplicity, we neglect the contribution of radiation. For cosmology on the brane (BC), Eq. (pM)) 
can be rewritten as 

D L° = 4-i 1 + z ) \ { [Vm(l + zf + n D p(z)/ Po ] [1 + 5(n m (l + zf + Q D p(z)/p )] y 1/2 (1 + 5f/ 2 dz , (25) 

-"0 Jo 

where the convenient parameter S — (po + p m )/2X has been introduced. 

IV. DARK ENERGY MODELS ON THE BRANE WITHOUT FINAL SINGULARITIES 

The following are examples of exact solutions corresponding to the Little Rip case, asymptotic de Sitter regime, 
and asymptotic breakdown. It will be assumed, for simplicity, that the universe consists of dark energy only. 

A. Little Rip case 

The simplest choice of EoS is to set f(rj) = a 2 — const. Then, 

* 1 \ ( v /2A 77 
3 -,cosh77j, %) = Vy^ 



0(77) = a exp ( — cosh 77 , £(77) = d — — ^ . (26) 



that is 

a{t) = a exp cosh | v/^" 2 * 



(27) 
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The dependence of the phantom energy density on time reads 

21 i^a 2 



p(<) = 2Asinh^ h/^Y* 



(28) 



therefore, a Little Rip occurs on the brane exponentially faster than in the usual FRW cosmology (see pjlQil)- The 
scalar field grows linearly with time, as 



4> = 0o + at. 

The potential of the scalar field is exponential, as expected in a string theory: 



V{(j>) = 2Asinh 2 



3 a( 



)) 



a 
Y 



(29) 



(30) 



The acceleration of the universe leads to an inertial force on any mass, m, as seen by a gravitational source separated 
by a comoving distance I, of the form 



jr. 

1 7,7 



il— = ml 



P , f(p) , P 2 , }{P)P 
6A 



3 2 6A 2A 
Following (18l |. it is convenient to define the dimensionless parameter 

a _ 2p + 3f(p)+p 2 /X + 3f(p) P /X 



F m = 6- 



apo 



(31) 



(32) 



where po is the dark energy density at present. 

The first two terms in Eq. ([31]) coincide with the inertial force in the FRW cosmology. From ([31]) one can conclude 
that the dissolution of bound structures in a FRW brane universe occurs essentially faster than in the usual FRW 
cosmology. The inertial force is a nonlinear function of the energy-density, which in our model reads 



p = 2X sinh 7* + Arsh 



Via 2 



2V2A 



For t >> 7 , one can simply take 



A 



P ~ g cx p(t£) 

and, therefore, the dimensionless inertial force grows exponentially for large times, namely 

A 



F 



^po 



exp(27i). 



(33) 



(34) 



(35) 



Performing similar calculations for the classical cosmology, one finds that the inertial force increases with time as 
t 2 , indeed 



- 3a 2 , ,n 

F in « (at) 2 . 

Po 



(36) 



These models can be matched with the latest data from the Supernova Cosmological Project. For illustration, let us 
rewrite the energy density as a function of the redshift. We have 



P(z) = Po 



l--h(l + Z ) 
Po 



(37) 



and, thus, if the value of the dimensionless parameter A = 3a 2 /on << 1, the energy density is, with good accuracy, 
a constant: our model actually mimics a cosmological constant [a, @. Therefore, it is trivially, in a way, that our 
dark energy model describes, for A << 1 with fi^ = 0.72 and f2 m = 0.28, the standard FRW cosmology and shares 
its celebrated matching with the data coming from astronomical observations . For the brane model one should also 



7 



define the parameter S in Eq. (|25|) . Consequently, at least in the domain S « 1, the model (|37|) on the brane is 
indistinguishable from ACDM cosmology. 

All structures disintegrate when the inertial force (|3"Tj) . dominated by dark energy, becomes equal to the forces 
which keep them bounded. For the system Sun-Earth, for instance, this happens when Fi n reaches a value ~ 10 23 . 
Let us choose, for an example, A = 0.01 and S = 0.01. In this case, the parameter a = (f^A) 1 / 2 /^ ~ 0.00624 
Gyr -1 . The time required for the disintegration of our solar system is 

t dls « 5 x 10 14 Gyr 

in FRW cosmology, but only of 

t dls a 9.4 x 10 5 Gyr 

in the brane model. The disintegration time slowly increases with the growth of the brane tension; to wit, for 5 = 10 -19 
the disintegration time becomes closer to the corresponding time in the FRW cosmology, being then tdis ~ 0.6 x 10 14 
Gyr. 

Summing up, we have presented above a realistic Little Rip cosmology on the brane coming from a scalar dark 
energy with an exponential potential, and have proven that the Little Rip regime sets up remarkably faster here than 
in the convenient FRW cosmology, owing to the presence of the brane tension. 



B. Asymptotic de Sitter regime 

We will now consider the situation corresponding to the case of a brane with negative tension. Let us choose 

f(rj) — a 2 cos 2 (77), a — const, (38) 

which translates into the following EoS 

p(p) = -p-a 2 (l+p/\) 2 . (39) 
Time, as a function of the parameter 77, reads 

'2|A^ 1/2 



a (tan 77 — tan 770) 



(40) 



and we see that, when 77 — > 7r/2, then t — > 00 while a — > const. That is, the universe expands in a quasi-de Sitter 
regime. The scale factor increases with time as 



a(t) = ao exp 



/3 2 



1 /2 

(l + (tan % + at/P) 2 ) -cos _1 77 



(41) 



In Eq. (j4Tj) the dimensionless parameter f3 2 — 2|A|/3a 2 has been introduced. When t — > 00, the scale factor grows 
exponentially with time, 

a(t) - exp [(|A|/6) 1/2 t 

Therefore, the effective cosmological constant is A e /j = |A|/2. The potential of the scalar field develops a complicated 
shape. The time dependence of the scalar field is found to be 



,.. n 1 — tan<I> 1 

<p(t) — — pln^ — ; - — — , < P = — arctan(tan ?7o exp{at/ p)), 



(42) 



where integration constants are omitted. The potential 



V($) = — cos 2 2$ + A(cos 2$ - 1) 



(43) 



tends to |A| with growing time. Note that, in the brane case, the maximal value of the inertial force (for A > 0) can 
considerably exceed the inertial force in the FRW cosmology, for the chosen EoS for dark energy. To illustrate this 
possibility, we consider a very simple model, with the following EoS 



f(p) = a 2 1 - 



Pf 



(44) 
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Observe that, for po < Pfi this equation describes phantom energy while, for po > pf, it corresponds to quintessence. 
The energy density-redshift relation reads 

Po 



p(z)=p f [1-(1- S ) (l + z) A5 °], 5 o = ^, A = 3a'/p . (45) 
When t — > oo the universe expands in a quasi-de Sitter regime. The scale factor grows exponentially with time, as 

a(t)~exp(y^V (46) 

The effective cosmological constant is then A e // = pf(l + S/Sq), and one concludes that A e // >> p/, if pf >> A. 
Also in this case pf >> A, for the dimensionless inertial force at t — > oo, one gets 



F 



2 Pf Pf >:> 

po 2A po 



the last expression being the inertial force for the FRW cosmology. 

We assume, for simplicity, that pf >> po >> a 2 . Using Taylor expansion, one can rewrite Eq. (|4"5|) as 

p(z) Aln(l + *)), 



(47) 



(48) 



with the same dimensionless parameter A as in Eq. ([571) . Therefore, as in the former case of the model with a 
Little Rip, the new model with EoS given by (|44)) , for 6 << 1, is also in perfect agreement with the most accurate 
observational data. Indeed, if we choose Po/pf = 10~ 13 and let 6 — 0.001, the value of the dimensionless inertial force 
on the brane is of the order of 10 23 (and thus, sufficient for the disintegration of our solar system) while it is only of 
some 10 13 in the case of ordinary FRW cosmology. 

As a consequence, on the brane a realistic mild phantom scenario (called a pseudo-Rip) is easy to realize. Owing 
to the brane tension, the effective "cosmological constant" in brane cosmology could be sufficiently larger than for 
FRW cosmology and, therefore, disintegration of bound structures may take place in the first case even if the second 
one remains safe yet. 



C. Asymptotic breakdown 



Asymptotic breakdown can be realized only for a brane model with negative tension and provided the extra 
dimension is time- like. If 77 — >■ 7r, the Hubble parameter tends to zero i.e., the universe bounces when the dark energy 
density has reached a sufficiently large value. Therefore, the expansion rate of the universe decreases with time. There 
are two possibilities: (i) the scale factor grows, albeit more and more slowly with time, or (ii) the scale factor tends 
to a constant value. For the first case, consider 



f(v) = ° 2 sin 2 7?, 



where the EoS in explicit form can be written as 



2a 2 
|A| 



P+^P 



(49) 



(50) 



Then, 



p / co t no\i 3 '/ 2 

Kv) = -(cot 77 -cot r/o ), a (77) = a — -|- 



and. therefore, 



a(t) = a 



+ (i+x 2 ) i ' 2 Y 2 ' 2 



cot • 



x = at/ (3 — cot 770. 



(51) 



The asymptotic behavior of this expression is power- like when t — > 00 (x — > 00), namely 

a(t) ~ t^l 2 . 
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If one chooses f(rj) = a 2 sin 2 T], then the universe expansion will asymptotically stop. After a rather simple calculation, 
we obtain for the scale factor 



a(t) = ao exp (3 2 arctan ^tan ^ exp(af / /3)j — /3 



2V0 
2 



Therefore, a — > aoexp[/3 2 (-7r — rjo)/2] when t — > 00. It is interesting to note that, at some point, the accelerated 
expansion of the universe turns into deceleration. As one can easily see from the following expression for the second 
derivative of the scale factor, 

a 9a sin 77 . 

- = — T\2 ( sm V + cos r ?) ' ( 52 ) 
a 4A Z 

this instant corresponds to 77 = 3ir/4. 

V. DARK ENERGY MODELS ON THE BRANE WITH FUTURE "SOFT" SINGULARITIES 

In this section we study brane dark energy models which lead to Type II or Type III future singularities and 
compare the corresponding predictions with the ones of the standard FRW cosmology [l9| . 

A. Big Freeze singularity (BFS, Type III) 

Let us consider this case for A > 0. Take 

f{rj) = a 2 cosh 2 rj (53) 
which leads, in explicit form, to the following EoS 

p-a 2 ^ + l)\ (54) 



V 



a, choice which yields a solution exhibiting two BFS: one in the past (rj — > —00) and another in the future (r) — > +00). 
Time is given by 

, . /2Atanh?7 

'W = Vt-^' (55) 

and the universe begins at U n = —-^ \J~^- and ends its existence at tf = ^ The expression for the scale factor 

eads 



a(t) = a,f exp 



A ( 3a 4 t 2 \ 1/2 ' 



3a 2 V 2A 



(56) 



where a/ is the final value of the scale factor. This dependence shows that the universe will contract in the interval 
t in < t < and then it will expand. 

Can a model like this be matched with actual observational data? Remarkably, the answer is positive. Indeed, after 
a simple calculation we obtain, for the dark energy density as a function of the redshift, z, 

l-(l/2 + n D J)Aln(l + s) 2 
P = P ° l + (l + 2^)A^ln(l + z)' A = 3 "/^ ( 5? ) 

Once more, we see that when 8 « 1 and A << 1, our model precisely mimics ACDM-cosmology, specifically in the 
range < z < 1.5. The time at the final singularity is 

In particular, if A = 0.1 and 8 = 0.01, then the time elapsed before the singularity is formed is of some 2 x 10 3 Gyr. 
Other models with similar properties can be readily constructed in an analogous way. 
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B. Sudden future singularity (SFS, Type II) 

Let us now discuss realistic models of dark energy which contain a Type II future singularity for the brane with a 
negative tension. To this end, we choose the following EoS 

f( r )) = —2T °<V*<k/2. (59) 

cos^(?7 — 77*) 

where j3 is a constant, either negative or positive. This EoS describes both phantom energy (when (3 > 0, < 770 < t/2) 
and quintessence (when j3 < 0, ir/2 < rj Q < n). The explicit function f(p) reads 

f{p)= (l-p/X*)* ' ^=2|A|sin 2 (J + ^)>|A|. (60) 

When p — > A*, the pressure diverges and a SFS is formed. The Hubble parameter remains finite but the acceleration 
of the universe diverges, a — > ±00. 

Let's consider the case when 77* is close to 0. With good accuracy one can derive the following relation 



p(z)n\\\{l- 



1/3' 

W\) + W\ {z+ ] 



(61) 



If (3 << po << |A| (phantom energy), this model mimics ACDM cosmology in the range of z corresponding to the 
SNe observations: 

p(z)&p (l- Aln(z + 1)), A = 3/3/p . 

The time remaining until the final singularity is formed can be calculated from Eq. (|20|). The instant the SFS occurs 
corresponds to 77 = it/2 + 77* w tt/2. Therefore, one has 

^ /2|A|\ 1/2 1 /tt % sin 2770 ^ 



and taking into account that |<5| << 1, and correspondingly 770 << tt/2, 

tsFS • 



2|A|\ 1/2 k il 1 



3 J (3 4 4 H \S\V 2 An D ' 

In particular, if 5 = —0.01 and A = 0.1, the remaining time is 1.5 x 10 3 Gyr. Acceleration of the universe turns into 
deceleration due to the negative tension of the brane. The universe acceleration as a function of the density is (we 
neglect here matter contribution and assume that dark energy dominates at present) 

a = -A l+ 2\) + 3(1- P /A*) a U + 2aJ- (63) 

When dark energy density reaches |A| the second term in this equation vanishes and then changes sign from "+" to 
"— ". Near the sudden future singularity (f(p) — > ±00), the expression for the inertial force (|3"Tj) reduces to 

p m±mp£. (64) 

Po 

In our case, p — > A*, Eq. (|ul|) yields 



and hence, due to the negative tension on the brane, we have a so-called "Big Crush" : at p —> A* > | A| the acceleration 
of the universe a — > — 00, instead of d — > 00 as for the usual FRW universe. 
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To finish, we can consider the scalar description of the model under investigation. For the scalar field, as a function 
of 77, we have 

/ 2|^| \ 1 / 2 

<j>{rj) = cf> } sin(»7 - rf), ^=1 — 1 , 
where the integration constant has been omitted. From here, the potential as a function of the scalar field is given by 

m = § — + m " ^ ,4>2f)1 ' 2 l) - (65) 

When (j> is small (77 — > in the phantom case, and 77 — > it for quintessence) this potential corresponds to a massive 
scalar held, namely 

V(<f>) 0«O. 
2<t> f 

Then, the sudden future singularity is here the singularity of the potential at <fi = <pf. In exactly the same way, 
comparison of other brane dark energy models, leading to soft singularities in the future with the predictions coming 
from the ordinary FRW dark energy, can be done. 



VI. CONCLUSION 



Several dark energy models on the brane, with various different scenarios of evolution have been presented. We must 
emphasize the remarkable property, shared by all these models, that they are able to mimic, under quite reasonable 
conditions, the standard ACDM cosmology at the present epoch but exhibit, otherwise, very different behaviors con- 
cerning their future evolutions, including asymptotic de-Sitter evolution, Little Rip behavior with eventual dissolution 
of bound structures, and various different future singularities (as the Big Rip one, Type II and Type III singularities). 

The future evolution of the universe is actually determined by the selected EoS for the dark energy component. The 
presence of a brane tension has been proven to lead to very interesting effects. In particular, owing to its influence 
the Little Rip occurs faster on the brane than in the usual FRW universe. Also, at a pseudo-Rip the effective value of 
the vacuum energy does not coincide with the dark energy density and, consequently, the disintegration of all bound 
structures can actually occur. Furthermore, a negative brane tension changes acceleration into deceleration and leads 
to a "Big Crush" of the universe, its existence ending, in this case, in a sudden future singularity. 

These results can be eventually used to discriminate the theoretical dark energy models when confronting them with 
precise cosmological data. The consequences of all these calculations are to be taken into account quite seriously, since 
they could actually compromise the degree of predictiveness commonly attributed to the more standard cosmological 
scenarios. 
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